Astronomy & Astrophysics manuscript no. capalbi 


February 5, 2008 


(DOI: will be inserted by hand later) 





Long-term monitoring of the X-ray afterglow of GRB 050408 with 

Swift/XRJ 

M. Capalbi 1 , D. Malesani 2 , M. Perri 1 , P. Giommi 1 ' 3 , S. Covino 4 , G. Cusumano 5 , V. Mangano 5 , T. Mineo 5 , S. 
O ■ Campana 4 , G. Chincarini 4 ' 6 , V. La Parola 5 , A. Moretti 4 , P. Romano 4 , G. Tagliaferri 4 , L. Angelini 7 - 8 , P. Boyd 7 , D.N. 
Burrows 9 , O. Godet 10 , J.E. Hill 7 ' 11 , J.A. Kennea 9 , F. Marshall 7 , PT. O'Brien 10 , and N. Gehrels 7 

O 



(N 



> 
in 

oo 

o 

O 

6 



13 



ASI Science Data Center, via G. Galilei, 1-00044 Frascati, Italy 
International School for Advanced Studies (SISSA-ISAS), via Beirut 2-4, 1-34014 Trieste, Italy 
ASI, Unita Osservazione dell'Universo, viale Liegi 26, 1-00198 Roma, Italy 
INAF - Osservatorio Astronomico di Brera, via E. Bianchi 46, 1-23807 Merate, Italy 

INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica, Sezione di Palermo, via U. La Malfa 153, 1-90146 Palermo, Italy 
Universita degli Studi di Milano-Bicocca, Dipartimento di Fisica, piazza delle Scienze 3, 1-20126 Milano, Italy 
NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA 

Department of Physics and Astronomy, The Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 21218, 
USA 

Department of Astronomy and Astrophysics, Pennsylvania State University, University Park, PA 16802, USA 
Department of Physics and Astronomy, University of Leicester, Leicester LEI 7RH, UK 

Universities Space Research Association, 10211 Wincopin Circle, Suite 500, Columbia, MD, 21044-3432, USA 



Received / Accepted 
Abstract. 

We present observations of the X-ray afterglow of GRB 050408, a gamma-ray burst discovered by HETE-II. Swift began ob- 
serving the field 42 min after the burst, performing follow-up over a period of 38 d (thus spanning three decades in time). The 
X-ray light curve showed a steepening with time, similar to many other afterglows. However, the steepening was unusually 
smooth, over the duration of the XRT observation, with no clear break time. The early decay was too flat to be described in 
terms of standard models. We therefore explore alternative explanations, such as the presence of a structured afterglow or of 
long-lasting energy injection into the fireball from the central GRB engine. The lack of a sharp break puts constraints on these 
two models. In the former case, it may indicate that the angular energy profile of the jet was not a simple power law, while 
in the second model it implies that injection did not stop abruptly. The late decay may be due either to a standard afterglow 
(that is, with no energy injection), or to a jetted outflow still being refreshed. A significant amount of absorption was present in 
the X-ray spectrum, corresponding to a rest-frame Hydrogen column density N H = 1.2^ x 10 22 cirr 2 , indicative of a dense 
environment. 
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1. Introduction 

T he Swift Gamma-Ray Burst Explorer dGehrels et alJ 
I2004I) . designed to study gamma-ray bursts (GRBs), has 
unique character istics all owing th e prompt observation o f GRB 
afterglows (see iBarthermv et all 120051 burrows etaD [2005; 
iRoming et all2005[ for a description of the Swift instruments). 
This has opened up the opportunity to study the previously un- 
explored early phases of their evolution. 

On the other hand, being fully devoted to GRB studies, 
Swift also has the capability to perform detailed, long-term 
monitoring of afterglows. To date, a number of GRBs have 
been observed up to ~ 1 month after the explosion (e.g. 
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XRF050416A; iMangano etaH Eo06). allowing a systematic 
study of the late phases of their evolution. Swifts rapid slew- 
ing capabilities also allows the prompt follow-up to alerts com- 
ing from other satellites, such as HETE-II and INTEGRAL. 
GRB 050408 was the first burst observed by Swift which 
was not autonomously-triggered by the Burst Alert Telescope 
(BAT) barthelmy onboard. 

GRB 050408 was detected by HETE-II at 16:22:50.93 UT 
on 2005 April 8 and localized by the Soft X-ray Camera (SXC) 
at RA(J2000) = 12 h 02 m 15 s , De c(J2000) = 10°52'01" (80" 
error radius, 90% containment: ISakamoto et ail 1200 5). The 
burst duration was Tgo ~ 34 s in the 7-40 and 7-80 keV 
bands, and Tgo ~ 15 s in the 30-400 keV band. Spectral 
analysis of the prompt emission showed that the fluence was 
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~ 1.4 x 10~ 6 erg cirT 2 and ~ 1.9 x 1(T 6 erg cm -2 in the 2- 
30 keV and 30^+00 keV energy bands, resp ectively. Using the 
classification defined bv lLamb etal] J200^ . GRB 050408 can 
there fore be classified as an "X-ray rich" GRB dSakamoto et alJ 
2005). 

On the basis of the HETE-II position, reported soon after 
the discovery via a GCN notice, a target of opportunity (TOO) 
was uploaded to Swift and the NFIs were pointed at the target 
about 42 min after the trigger. The Swift XRT detect ed an un- 
catalogued, fading source inside the SXC error box dwells et alJ 
120051). which was later refined with an error circle of 5" radius 
(Chincar ini et alJ20 05). Following the initial observation, XRT 
continued monitoring the afterglow for several weeks, leading 
to a long, well-sampled light curve. 

Several follow-up observations were performed at other 
wavelengths. An optical counterp art was discov ered inside the 
SXC error circle dde Ugarte Postigo et alJl2005l . using the lm 
and 6m telescopes at the Special Astrophysical Observatory. 
This obje ct was later se en to fade, confirming its afterglow na- 
ture ( Huang et al. 2005 ). The astrometric position was provided 
bv lChen et al J ( 120051) . who located the source at the coordinates 
RA(J2000) = 12 h 02 m 17.328, Dec(J2000) = +10°51'09'/47 
(0725 error radius). A spectrum of the optical afterglow ob- 
tained with the LDSS-3 instrument on the Magellan/Clay tele- 
scope showed the presen ce of emission and ab sorption lines 
at a redshift z = 1.236 dBerger et al.ll2005albl) . The Gemini 
telescope found a consistent value of z = 1.2357 + 0.0002 
dProchaska et alJl2005_ iFolev et alJl2006l) . Many other optical 
telescopes, including the 3.6m Telescopio Nazionale Galileo 
(TNG), detected and monitored the optical afterglow. The anal- 
ysis of the TNG data, taken simultaneously with the XRT 
observation, is presented in a separate paper dCovino et alJ 
120061). Some of the a fterglow properties have been discussed 
bv lFolev et alJ (2006). At the location of the afterglow, UVOT 
detecte d a faint, low-signific ance source in the coadded f/-band 
image ( Holla nd" et alJl2005l) . No detection was possible in the 
other optical and ultraviolet filters. Very Large Array radio ob- 
servations at 8.5 GHz revealed no sources at the afterglow po- 
sitio n, down to a 2-cr upper limit of 74 yidy on 2005 Apr 9.26 
TJT dSoderbergl2005l) . 

In the following sections we report a detailed analysis of 
the XRT follow-up observations. In Sect. |2] we describe the 
XRT data reduction and analysis, and in Sect. [5] we discuss 
the results. A summary of our work is reported in Sect. |4] 
Throughout the work, we follow the notation F v (t, v) oc r"v^ 
for the time and spectral dependency of the flux. The times 
quoted are with respect to the HETE-II trigger time. 

2. XRT data analysis 

XRT observations of the GRB 050408 field started on 2005 
April 8 at 17:05:17 UT (2545 s after the HETE-II trigger); 
data collection in photon counting (P C) mode s t arted a few 
seconds later f see .Hill et ail (120041) and Iffill et al.ld2005l) for a 
description of XRT readout modes). Swift subsequently contin- 
ued monitoring the GRB field at later times, collecting a total 
of 12 observations over a period of 38 d. The total exposure 
time in PC mode was ~ 223 ks. Occasionally, due to the high 



Table 1. X-ray afterglow light curve. The first column reports 
the center of the time bin, expressed in seconds since trigger. 



Time 


Bin size 


Count rate 


(s) 


(s) 


(10~ 2 count s _1 ) 


2612 


120 


23+5 


2732 


120 


30±5 


2852 


120 


29±5 


6302 


500 


19±3 


6802 


500 


15±2 


7302 


500 


16.6+1.8 


7802 


500 


13.0±2.0 


8302 


500 


14.8+1.7 


13302 


500 


9.7+1.5 


13802 


500 


10.8+1.5 


14302 


500 


8.8±1.4 


19202 


900 


6.1±1.2 


20102 


900 


6.7±1.0 


25052 


1000 


6.0+1.2 


26052 


1000 


5.0+1.0 


30482 


2660 


3.5+0.7 


40002 


10700 


2.7+0.6 


50702 


10700 


3.1+0.6 


61402 


10700 


2.7 ±0.5 


281121 


52356 


0.53±0.12 


418938 


80974 


0.47+0.13 


618410 


100110 


0.24+0.05 


1492073 


353910 


0.044±0.010 


1842640 


173812 


0.050±0.010 


2476438 


225974 


0.035±0.012 



background, the XRT switched into windowed timing mode. 
However, the source count rate was very low in these frames, 
and therefore these data were not included in our analysis. 

The XRT data from the SDC archive at the NASA Goddard 
Space Flight Center were processed with the XRTDAS 1 software 
package (v. 1.7.1) developed at the ASI Science Data Center. 
Calibrated and cleaned level-2 event files were produced with 
the xrtpipeline task, applying the standard screening crite- 
ria: frames with a CCD temperature greater than -47 °C were 
rejected, and bad pixels and bad aspect time intervals were 
eliminated. 

2.1. Image analysis 

The 0.3-10 keV PC mode image of the field was analyzed 
with the XIMAGE package (v. 4.3). A previously uncatalogued, 
bright X-ray source was clearly detected inside the HETE-II 
error circle, at the coordinates RA(J2000) = 12 h 02 m 17 s .29, 
Dec(J2000) = 10°51'11'.'4 (3'.'5 error radius; 90% contain- 
ment). This position and its error were evaluated taking into ac- 
count the latest XRT boresight calibration dMoretti et alJ20 06). 
The XRT coordinates are 60" and 1795 away from the HETE-II 
error box center and the optical afterglow position, respectively. 
Another non-fading source with count rate (2.9 + 0.1) x 10" 3 
count s _! is present inside the HETE-II error circle, at a dis- 
tance of « 38" from the afterglow. To avoid contamination 
from this source, the events for the temporal and spectral analy- 

http : //swift . gsfc.nasa.gov/docs/swift/analysis/xrt_swguide_vl_2 .pdf 
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Table 2. Best-fit parameters with their associated 68% confidence intervals, for different functional forms: single power law (PL), 
sharply broken power law (2PL), smoothly broken power law (SBPL), and Beuermann fit ("B : iBeuermann et afll 1 999l) for two 
values of the smoothness parameter s. 



Model 


)G (dof) 




on 




a<2 






fb(d) 






Best 


68% C.I. 


Best 


68% C.I. 




Best 


68% C.I. 


PL 


2.50 (23) 


0.99 


0.97 < a, < 1.01 












2PL 


1.0 (21) 


0.80 


0.49 < o-i < 0.67,0.76 < or, < 0.84 


1.40 


1.25 < a 2 


< 1.55 


2.2 


0.12 <t b < 0.18, 1.3 <t b < 3.0 


SBPL 


0.92 (21) 


0.46 


-0.12 <a t < 0.78 


1.24 


1.11 < a 2 


< 1.91 


0.21 


0.05 <t b < 6.6 


B, s = 0.1 


0.89 (21) 


0.21 


a, < 0.31 


2.05 


1.81 < a 2 


< 2.89 


unconstrained 


B, s = 0.5 


0.89 (21) 


0.13 


ai < 0.75 


1.51 


1.15 < a 2 


< 3.10 


unconstrained 



10 -io 




10-' 10° 10' 

Time since burst (d) 

Fig. 1. The 0.3-10 keV X-ray light curve of the afterglow 
of GRB 050408 (PC mode data only). The solid and dashed 
lines show the best-fit models assuming a sharply- (2PL) and 
smoothly- (SBPL) broken power law, respectively. 



sis were selected from a circular region of 8 pixels (19") radius 

centered at the afterglow position, which contains about 77% 

I i I 

of the photons at 1.5 keV (Moretti et al. 2005). The appropriate 

ancillary response file was used to correct for the PSF losses. 



2.2. Temporal analysis 

For the temporal analysis, the standard grade selection for 
PC mode (grades from to 12) was adopted, in order to maxi- 
mize the statistics. Only the 0.3-10 keV energy range was con- 
sidered. The background in the extraction region was evaluated 
by comparing several source-free boxes in the field and, since 
it was stable during the whole observation, a constant level 
of 1.3 x 10~ 4 count s _1 was subtracted from the light curve. 
Moreover, the contribution of the nearby serendipitous source 
was evaluated computing the number of photons falling inside 
the afterglow extraction region and the afterglow count rate was 
corrected for this contamination. The light curve was binned to 
ensure a minimum of 20 counts per bin. The data are reported 
in Table ^ The count rate was converted to an unabsorbed 0.3- 



10 keV flux using a conversion factor of 1.0 x 10 erg cm"" 
count -1 , obtained using the results from the spectral analysis 
(Sect 2.3). The light curve is shown in Fig.^points). A gap of 
approximately two days is present in the data due to the close 
proximity of the Moon preventing the NFIs from observing. 

We initially fit the afterglow decay with a single power 
law, obtaining an unacceptable^ = 2.5 (23 degrees of free- 
dom, d.o.f.). Looking at the residuals, a systematic deviation 
of the data from the model was apparent, showing a continu- 
ous steepening of the light curve with time. We thus fitted the 
data with a broken power law (F oc r"' for t < fb, F oc r" 2 
for t > fb), where ct\ and ai are the early- and late-time de- 
cay slopes, respectively, and fb is the break time. The best fit 
(Fig-E solid line) significantly improved, yielding^ = 1.0 
for 21 d.o.f. (null hypothesis probability of 6.6 x 10~ 5 , accord- 
ing to an F -test). The best fit parameters, together with their 
confidence intervals, are reported in Table|5] Errors were com- 
puted leaving all parameters free to vary. The fit, and in partic- 
ular fb, is not well constrained. Fig.|2](upper panel) shows the 
behaviour of as a function of fb. Two minima are apparent 
at approximately the same level (with the larger va lue of h be- 
ing sl ightly preferred). In an independent analysis. iFolev et all 
(2006) indicated the lower value as the best fit. The reason 
for the discrepancy is possibly due to minor differences in 
the reduction process (e.g. background subtraction, data bin- 
ni ng). However, ou r resul ts ar e consistent w i th tho se derived 
bv IChincarini et al J i2QQ(h and lNouseket al.1 i2QQ(ii (who an- 
alyzed only the first part of the light curve). The smoothness 
of the light curve prompted us to attempt different functional 
forms, using the so-called smoothly-broken power law (SBPL) 
model: F(t) oc l/[(f/f b ) Q ' 1 + (f/fb)" 2 ]- This provided a slightly 
better fit Cr r 2 = 0.92 for 21 d.o.f.; dashed line in Fig.Q), but 
fb could not be constrained (Fig. |2 middle panel). Similarly, 
fb could not be constrained with a fit to the Beuermann model 
(F(t) oc [(f/f b ) OT > + (f/f^r'^ fBeuermann etakll 19991) . The 
parameter s (the smoothness index) controls how fast the tran- 
sition between the two phases happens: large values of s imply 
a sharp break. Leaving all parameters free, small values of s 
(i.e., a smooth break) were systematically preferred (s = 0.1- 
0.5) but the fit did not converge. Fixing s to a set of definite 
values allowed the confidence intervals for the parameters to 
be determined (Fig. [2 lower panel). 

Finally, in consideration of the fact that the values of the 
last three data points of the light curve are consistent within the 
errors, we checked for the presence of a persistent X-ray emis- 



4 



M. Capalbi et al.: Swift/XKl observations of the X-ray afterglow of GRB 050408 



1.5 
1.0 h 



1.1 

1.0 E- 
0.9 'r 



0.8 '- 
lO" 2 



1 II I I I I I I I 1 1 1 I I I I I I 1 1 — I I I I I I I 



2PL 

i i i 1 1 1 1 1 1 1 i i 1 1 m l 1 — i i 1 1 n il 




i i i . i i n i| 1 i i 1 1 iii| 1 i i 1 1 1 1 1 , | h 
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s=0.1 



s=0.5 
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10-' 10° 
Break time (d) 



10' 



Fig. 2. Reduced x 1 as a function of the break time fb for 
three functional forms: from top to bottom, joint broken power 
law (2PL), smoothly broken power law (SBPL), and the 
Beuer mann fit, for different values of the smoothness param- 
eter s JBeuermann et alll9 99). 




Energy (keV) 

Fig. 3. The average 0.3-10 keV X-ray spectrum of the 
GRB 050408 afterglow (top panel), together with the best-fit 
absorbed power-law model (with Galactic and host absorption). 
The bottom panel shows the ratio between the data and the best- 
fit model. 



sion coming from the afterglow host. We found that adding a 
constant component to the broken power law model (2PL) does 
not yield a significant improvement to the fit (F-test probability 
of a chance improvement ~ 80%). 



2.3. Spectral analysis 

Events were extracted for the spectral analysis from the 
same circular region used for the temporal analysis, but a more 
strict selection on event grades was applied (grades 0-4, i.e. 



single- and double-pixel events), in order to achieve better 
spectral resolution. The spectrum was binned to ensure a mini- 
mum of 20 counts per energy bin. Channels below 0.3 keV and 
above 10.0 keV were excluded. We first produced an average 
spectrum obtained by summing all the available observations. 
A spectral fit was done using the XSPEC package (v. 1 1.3.1). We 
modelled the spectrum with an absorbed power law of spec- 
tral index /3 and Hydrogen column density Nu (assumed to be 
at redshift z = and with the heavier elements fixed at so- 
lar abundances). The best-fit parameters are shown in Table[5] 
The value of Nh was found to be significantly high er than the 
Galac tic value (N H = 1.74 x 10 2Q cm" 2 . iDickev & Lockmar] 
Il990t N H = 1.54 x 10 20 cirr 2 . lKalberla et al.l2005l) . We fhere- 
fore included an additional absorption component at the red- 
shift of GRB 050408 (model ZWABS in XSPEC), keeping values 
of the redshift (z = 1.2357) and the Galactic column density 
(N H = 1.74 x 10 20 cm" 2 at z = 0) frozen. The fit provided 
N H = 1-2!^ x 10 22 cm" 2 for the redshifted absorber. The 0.3- 
10 keV XRT spectrum is shown in Fig. [3] together with the 
best fit absorbed power law model. Our results are consistent 
with the analysis reported by Folev et al 1 d2006l) . The average 
unabsorbed 0.3-10 keV flux, calculated between ~ 2.5 ks and 
~ 67 ks after the trigger, is 7.1 x 10" 12 erg cm" 2 s" 1 . 

We looked for spectral variability across the observation. 
Due to the limited statistics at late times, we split the data into 
two bins, performing a power law fit (with absorption both in 
the Milky Way and in the host galaxy) to both sections. The 
best-fit parameters were consistent with those found for the av- 
erage spectrum (Table[3j. For the second part of the light curve, 
we had to freeze the column densities (which are not expected 
to change) to the values measured before the break. As can be 
seen from Table [5] the results from the two phases are consis- 
tent. To further test this result, we studied the time evolution 
of the hardness ratio, computed as the ratio between the count 
rate in the 2-10 and 0.3-2 keV bands. No significant variability 
was apparent. 

3. Discussion 

The Swift observations allowed a detailed study of the 
X-ray light curve of GRB 050408, extending for more than 
three decades in time. This in turn led to the detection of 
a very smooth break in the temporal decline, which would 
have been more difficult to detect with less coverage. The ini- 
tial portion of the light curve presents a shallow decay, fol- 
lowed by a steepening to a more conventional slope. This 
is consistent with the typical behaviour observed in sev- 
eral Swift GRBs dChincarinietal.lf2006HNouseketal.ll2006l 
lO'Brien et alJ2005h . assuming that the early rapid decline may 
have occurred during the 42 min between the trigger and the 
first obse rvation. The spectral slop e is also typical of GRB af- 
terglows (De Pasquale et al. 2005). 

What is less typical is the smoothness of the transition be- 
tween the fiat and steep phases. Despite the long-duration cov- 
erage, there is no clear evidence of a well-defined break time, 
fb. Indeed, fits to the light curve allow a very broad range 
for fb. Moreover, fits with a smooth transition are preferred. 
The behaviour has been observed in very few afterglows. Most 
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Table 3. Results of the spectral fits to the X-ray afterglow of GRB 050408. An asterisk indicates a frozen parameter. The redshifts 
for the Galactic and host absorption were kept fixed at z — and z = 1.2357, respectively. The errors are at 90% confidence level 
for one interesting parameter. 





Galactic absorption 


Galactic + host absorption 


Time range 




Nn 
10 21 cm- 2 


Xi (d.o.f.) 




Galactic Nh 
10 21 cm- 2 


Host N H 
10 21 cm- 2 


Xt (d-o.f.) 


All data 


1.3 ±0.2 


2.9 ±0.6 


1.4 (39) 


1.1 ±0.1 


0.174* 


l2 + _ 4 3 


1.5 (39) 


First part 


1.3 ±0.2 


3.0 ±0.7 


1.2 (32) 


1.2 ±0.2 


0.174* 


15 ±4 


1.3 (32) 


Second part 


1 2 +04 


3.0* 


1.1 (6) 


1.0 ±0.3 


0.174* 


15* 


0.9 (6) 



optical light curves, when fitted with the Beuer mann mo del, 
have smoothness para meters as large as s = 10 JCovino et alJ 
l2003UZeh et all2005h . In the X-ray bandpass, the statistics are 
usually worse, but joint power laws usually fit the data well. 
Without the long-term coverage provided by Swift, the curva- 
ture of this break would have easily been missed. 

The early part of the light curve (t < ft,) is too shallow to 
be explained in terms of the stan dard afterglow models (e.g. 
ISari et al.ll998HPanaitescu & Kumarl2000h . either considering 
a homogeneous (ISM) or a wind-shaped surrounding medium. 
The early-time slope is a\ < 0.84 for all models (see Tables- 
Given the spectral slope 2 f}\ = 1.2 ± 0.2, we have a\ -3/3\/2 < 
-0.96, which is inconsistent with the expected values: (ISM, 
v < v c ), +1/2 (wind, v < v c ), or -1/2 (v > v c ), where v c is 
the cooling frequency. Fast-cooling models are also ruled out, 
since a spectral index of fi\ = 0.5 would be expected. 

There are several ways to explain a shallow slope. For 
example, if the cooling of electrons is dominated by inverse 
Compton rather than synchrotron losses (that is, the Compton 
parameter is Y » 1), the decay above v c is expected to 
be flatter than in the sta ndard case (e.g. Sari&Esin 2001; 

In I 
Fanaitescu & Kumar 2000). In fact, the cooling frequency is 

v c = v^(l + Y)~ 2 (where vj! is the synchrotron cooling fre- 
quency), and it is easy to show that the flux above v c is F v (t) = 
Fl(t){\ +Y) , where F£ is the synchrotron flux in the case of no 
Compton losses. Since the Compton parameter decreases with 
time (at least in the slow cooling regime), when Y » 1 the 
decay law is flatter than in the case of pure synchrotron. The 
steepening of the light curve pinpoints the time at which Y — 1 
(so that F v ~ Fl). The smoothness of the transition would im- 
ply that the decrease of Y with time is very slow, possibly not 
following a power law. Given the soft X-ray spectral index, the 
Compton component should be confined to energies above the 
XRT range. 

Observations in the optica l show that this band likely lies 
below the cooling frequency JFolev et alJl2006r ICovino et alJ 
2006). The light curve at these frequencies may have a de- 
cay similar to that in the X-ray band, thus disfavoring this as- 
sumption. However, the poorly constrained X-ray decay and 
the presence of significant extinction in the optical makes more 
stringent comparisons difficult. 

An alternative way to explain the earl y flat decay is to 
have an outflow with an angular structure jRossi et afl r2002: 
IZhang & Meszarosll2002l) . that is a jet with an energetic core 

2 We denote with ji\ and fii the spectral slopes before and after the 
break, respectively. 



and less intense wings. As time elapses, the fireball Lorentz 
factor, F, decreases so that a larger portion of the jet becomes 
visible to the observer due to relativistic aberration. If the jet is 
viewed off-axis (at an angle # v ), the core contribution becomes 
more and more important, so that its increased emission partly 
compensates for the flux decline, leading to a slower decay. In 
this case, the light cu r ve bre ak would happen when F ~ 1/§ V . 
IPanaitescu & Kumarl J2003I) modeled the observed emission 
from off-axis structured jets, showing that flat {a\ « 0.5) early- 
time slopes are possible for & w /§ c ~ a few, where & c is the 
angular extension of the j et core. In this cas e, however, a sharp 
break is predicted (s > h lRossi etall2 004). Thus, the smooth- 
ness of the observed transition may indicate that the energy 
angular profile of the jet was not a pure power law as function 
of the off-axis angle. 

The final possibility is that the early shall ow decay is due 
to delayed energy injection in the fireball (e.g. | Dai & LvJl998l: 
IPanaitescu et aljfl998l IZhang & Meszarosll200ll) . This can be 
achieved in two ways: long-lasting energy emission from the 
central engine, or refreshed shocks from slow sh ells catc hing 
the leading fireball after it has decelerated (see e.g. lZhang et al. 
20061 for a general overview). Even if the latter model may 
be preferred on theoretical grounds (no extended activity is re- 
quired), the two scenarios are difficult to distinguish observa- 
tionally. Both can account for light curves as flat as F(t) oc t°. 
The smoothness of the break may indicate that, in this case, the 
additional energy injection did not follow a power law in time, 
but had a more complex history. 

However, an effective index q can be introduced by 
parametrizing the energy injection as £ oc r q . A viable so- 
lution is found with q < 0.70 (for a\ < 0.84), assuming v > v c 
(either ISM or wind). The case with V{ < v < v c (v; being 
the synchrotron injection frequency) is acceptable only in the 
ISM case (providing q < 0.48), while the wind case is excluded 
since q < would be required. Moreover, if v < v c , a very large 
p = 3.2+0.2 would be implied, so that we favor the case v > v c . 
Values of the or der q & 0.5 have been inferred also fo r other 
GRB afterglows (Zhan g et alJ2006HNousek et all2006l) . If en- 
ergy was supplied to the fireball through slow shells impacting 
the decelerated fireball, and parametrizing their Lorentz factor 
distribution as M(> y) oc y~ r (M(y) being the mass of shells 
with Lorentz factor larger than y), the above values for q cor- 
respond to r > 1.9 (ISM, v > v c ), r > 2.7 (wind, v > v c ), 
and r > 2.7 (ISM, v\ < v < v c ). In the context of the injec- 
tion model, two possibilities are viable to explain the observed 
break in the light curve. First, % may simply identify the end 
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of the energy injection process. The decay after would then 
correspond to a standard isotropic afterglow. Indeed, the ob- 
served ao ~ 1 .5 and fa » 1 .0 satisfy, within the errors, several 
closure relations (again, only the wind case with v, < v < v c 
is unfavored). An alternative explanation is that the injection 
did not stop at ft,, and the stee pening was due to a jet effect 
jRhoadsll999HSari et al.ll999t) . The data (Table are consis- 
tent with a late-time decay 0-2 < 2 (particularly if the break 
time was early), which cannot be accounted for in the standard 
model, where 0-2 — p t 2 is predicted. Energy injection con- 
tinuing after the jet break would naturally make the decay shal- 
lower. Indeed, there are examples where energy injection lasted 
for considerable durations (e.g. XRF 050406, where this phase 
lasted for > 10 6 s; Ro mano et al.l2006T) . A detailed analysis of 
the dynamics of a relativistic jet being refreshed is beyond the 
scope of this work, however we envisage this solution as quali- 
tatively possible. Finally we note that, within the uncertainties, 
our data do not exclude steeper slopes (a>i > 2), so that en- 
ergy injection is not strictly required after the break. However, 
it would be coincidental that the cessation of the refreshing and 
the jet break happened at the same time. 

The long-lasting injection episode implies th at the fire- 
ball energy increased significantly. iNousek et alJ d2006t) give 
several recipes to estimate the fractional energy increase, /. 
Limiting ourselves to the regime v > v c , and conservatively 
only taking the evolution where t < % into account, we have 
/ = (?b/ f start)*, where f start < 2500 s is the start of the injection 
phase, k = 2Aa/ ( 1 +/?) « 0.5 and Aa is the difference between 
the observed decay slope and the one which we would see in 
the case of no injection (a = 3/3/2 - 1 /2). With these numbers, 
and taking % ~ 10 5 s, we find / > 2 -f 25, so that the fractional 
energy increase was quite significant. In the case where v < v c 
(which we do not favor), / would be even larger. 

There is another remarkable feature concerning 
GRB 050408, namely its large rest-frame column density, 
of the order of 10 22 cirT 2 (computed assuming solar abun- 
dances). Such a value is typical of giant molecular clouds 
jReichart & Pried l2002h . Based on BeppoSAX and XMM- 
Newton data, many auth ors have reported column densities as 
large as 10 22 cm~ 2 (e.g.lO alama & Wiiersll200lHWatson et alJ 
2002J IStratta et alJl2004 be Luca et alibOOaT rhis was later 
confirmed by Campan a et alJ (|2006) using Swift-XKT data, 
showing that excess column density is a common feature 
among GRB afterglows. So, it is not surprising to find such a 
value. It is however more difficult to explain the relative bright- 
ness of the optical afte rglow. Using the Galactic gas-to-dust 
ratio JPredehl & Schmittll 19951) . the measured column density 
would correspond to Ay ~ 8 mag. This is not consistent 
with the detection of an optical counterpart (the observed 
R band would su ffer 14 mag of extinction at z = 1.2357). 
iFolev et all J2006). from optical spectroscopy and photometry, 
estimated Ay ~ 0.5-1 mag, assuming a SMC-like extinction 
curve. The small amount of optical extinction compared to 
the X-ray absorption was already noticed in previous cases 
( Galama & Wiienl200lllStratta et all2004 . and ratios Ay/Nn 
as low as ten times less than the Milky Way value have been 
reported. This implies either different dust optical properties 
or a low dust-to-metals ratio. 



In the former case, the estimation of the dust content from 
the afterglow spectral properties may be incorrect. For exam- 
ple, from the analysis of the absorbing element abundances in 
the GRB 020813 afterglow. ls"avaglio & Fall J2004 derived an 
amount of dust larger than inferred from the analysis of the 
continuum spectral sha pe, thus indi cating an anomalous ex- 
tinction curve. However, Fo ley et alJ (|2006) detected Titanium 
overabundance in the optical spectrum of the GRB 050408 af- 
terglow. Since this element is highly refractory, low dust con- 
tent (or a different dust composition) is inferred for this line of 
sight, independently on its transmission properties. This would 
therefore leave us with the second possibility, namely an in- 
trinsically low dust-to-metals ratio. Such composition may be 
a property of the GRB formation environments, where the in- 
tense UV radiation field from the young, hot stars may ham- 
per dust formation. Also, the young stellar age of GRB hosts 
may imply that dust has not yet had time to form \ Watson et alJ 
2006). Finally, a low dust content may be the direct effect of 
the burst explosion, which is able to subli mate dust grains up 
to ~ 100 pc from th e explosion site dWaxman & Drainel boOO: 
iFruchteret all200ll) . 

4. Conclusions 

GRB 050408, detected by HETE-II, was monitored by 
Swift-XKT over a period of 38 d. This allowed a detailed 
study of the X-ray light curve, extending over more than three 
decades in time. A very smooth break was apparent ~ 2 x 10 5 s 
after the trigger, separating a shallow decay phase from a 
steeper one. The transition was extremely gradual, with no 
definite break time. This is in contrast to the usual behaviour 
observed for X-ray and optical afterglows, where a sharp and 
abrupt break is usually observed. On the other hand, the spec- 
tral properties of this afterglow were not unusual, and no spec- 
tral variability was observed. The X-ray spectrum showed a 
large Hydrogen column density (Nh ~ 10 22 cm" 2 ), signifi- 
cantly in excess of the Galactic value. This may be a com- 
mon feature among GRB afterglows JCampana et all [2006; 
Galama & WiierslEoOlUStratta et alJl2004. po ssiblvindicat- 
ing dense environments (Reichart & Pricel2002l) . This large X- 
ray absorption was, however, acco mpanied by a relatively small 
optical extinction llFolev et all2 006: Covino et al. 2006). 

The first portion of the afterglow light curve is too shal- 
low to be explained in terms of standard afterglow models. 
Alternative solutions were considered, including the possibil- 
ity of a structured jet observed off-axis, or the presence of 
Compton radiation affecting the decay above the cooling fre- 
quency. Our data are also consistent with late injection of en- 
ergy from the central engine lasting for considerably longer 
than the GRB explosion. In this case, the break may pinpoint 
the end of the injection activity. Alternatively, the steepening 
may be a standard jet break. In this case the injection has to 
continue even during the steep phase. In both cases, the energy 
supplied by the central engine is considerable, being larger by 
a factor of > 2 -5- 25 than the initial amount. This may pose 
serious efficiency problems to GRB radiation models. 

Acknowledgements. We thank Don Lamb and Nat Butler for useful 
discussions, Ryan Foley and Dave Pooley for providing us their X-ray 



M. Capalbi et al.: Swift/XKl observations of the X-ray afterglow of GRB 050408 



7 



light curve. We thank F. Tamburelli and B. Saija for their work on the 
XRT data reduction software. We acknowledge financial support by 
the Italian Space Agency (ASI) through grant I/R/039/04 and through 
funding of the ASI Science Data Center. DM and OG gratefully ac- 
knowledge MIUR and PPARC funding. 



References 

Barthelmy, S.D., Barbier, L.M., Cummings, J.R., et al. 2005, 

Space Sci. Rev., 120, 143 
Berger, E., Gladders, M., & Oemler, G. 2005a, GCN Circ. 3201 
Berger, E., Kulkarni, S.R., Fox, D.B., et al. 2005b, ApJ. 634, 

501 

Beuermann, K., Hessman, F.V., Reinsch, K. et al. 1999, A&A, 
352, L26 

Burrows, D.N., Hill, J.E., Nousek, J.A., et al. 2005, 

Space Sci. Rev., 120, 165 
Campana, S., Romano, P., Covino, S., et al. 2006, A&A, 449, 

61 

Chen, H.-W., Green, P.J., Bloom, J.S., & Prochaska, J.X. 2005, 

GCN Circ. 3199 
Chincarini, G, Romano, P., Campana, S., et al. 2005, GCN 

Circ. 3209 

Chincarini, G., Moretti, A., Romano, P., et al. 2005, 

astro-ph/0506453 
Covino, S., Malesani, D., Tavecchio, E, et al. 2003, A&A, 404, 

L5 

Covino, S., Palazzi, E., Malesani, D., et al. 2006, A&A, sub- 
mitted 

Dai, Z.G, & Lu, T. 1998, A&A, 333, L87 
De Luca, A., Melandri, A., Caraveo, P., et al. 2005, A&A, 440, 
85 

De Pasquale, M., Piro, L., Gendre, B., et al. 2005, A&A, 455, 
813 

de Ugarte Postigo, A., Komarova, V., Fathkullin, T., et al. 2005, 

GCN Circ. 3192 
Dickey, J.M., & Lockman., F.J. 1990, ARA&A, 28, 215 
Foley, R., Perley, D.A., Pooley, D., et al. 2006, ApJ, 645, 450 
Fruchter, A.S., Krolik, J.H., & Rhoads, J.S. 2001, ApJ, 563, 

597 

Galama, T.J., & Wijers., R.A.M.J. 2001, ApJ, 549, L209 
Gehrels, N., Chincarini, G, Giommi, P., et al. 2004, ApJ, 61 1, 
1005 

Hill, J.E., Burrows, D.N., Nousek, J.A., et al. 2004, Proc. SPIE, 
5165, 175 

Hill, J.E., Angelini, L., Morris, D.C. et al. 2005, Proc. SPIE, 
5898, 325 

Holland S.T., Capalbi, M., Morgan, A., et al. 2005, GCN Circ. 
3227 

Huang, K.Y., Ip, W.H., Kinoshita, D., & Urata, Y. 2005, GCN 
Circ. 3196 

Kalberla, P.M.W., Burton, W.B., Hartmann, D., et al. 2005, 
A&A, 440, 775 

Lamb, D.Q., Donaghy, T.Q., & Graziani, C. 2005, ApJ, 620, 
355 

Mangano, V., La Parola, V., Cusumano, G, et al. 2006, ApJ, in 
press (astro-ph/0603738! 



Moretti, A., Campana, S., Mineo, T., et al. 2005, Proc. SPIE, 
5898, 360 

Moretti, A., Petri, M., Capalbi, M., et al. 2006, A&A, 448, L9 
Nousek, J., Kouveliotou, C, Grupe, D., et al. 2006, ApJ, 642, 
389 

O'Brien, P.T., Willingale, R., Osborne, J., et al. 2005, ApJ, 647, 
1213 

Panaitescu, A., Meszaros, P., & Rees, M.J. 1998, ApJ, 503, 314 
Panaitescu, A., & Kumar, P. 2000, ApJ, 543, 66 
Panaitescu, A., & Kumar, P. 2003, MNRAS, 592, 390 
Predehl, P., & Schmitt, J.H.M.M. 1995, A&A, 293, 889 
Prochaska, J.X., Bloom, J.S., Chen, H.-W., Foley, R.J., & Roth, 

K. 2005, GCN Circ. 3204 
Reichart, D.E., & Price, PA. 2002, ApJ, 565, 174 
Rhoads, J.E. 1999, ApJ, 525, 737 

Romano, P., Moretti, A., Banat, PL., et al. 2006, A&A, 450, 59 
Roming, P.W.A., Kennedy, T.E., Mason, K.O, et al. 2005, 

Space Sci. Rev., 120, 95 
Rossi, E.M., Lazzati, D., & Rees, M.J. 2002, MNRAS, 332, 

945 

Rossi, E.M., Lazzati, D., Salmonson, J., & Ghisellini, G. 2004, 

MNRAS, 354, 86 
Sakamoto, T., Ricker, G.R., Atteia, J.-L., et al. 2005, GCN Circ 

3189 

Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 

Sari, R., Piran, T., & Halpern, J.P 1999, ApJ, 519, L17 

Sari, R., & Esin, A.A. 2001, ApJ, 548, 787 

Savaglio, S., & Fall, S.M. 2004, ApJ, 614, 293 

Soderberg, A.M. 2005, GCN Circ. 3210 

Stratta, G, Fiore, F, Antonelli, L.A., Piro, L., & De Pasquale, 

M. 2004, ApJ, 608, 846 
Watson, D., Reeves, J.N., Osborne, J.P, Tedds, J. A., O'Brien, 

PL, Lomas, L., & Ehle, M. 2002, A&A, 395, L41 
Watson, D., Fynbo, J.P.U., Ledoux, C, et al. 2006, ApJ, in press 

(astro-ph/05 10368) 
Waxman, E., & Draine, B.T. 2000, ApJ, 537, 796 
Wells, A.A., Abbey, A.F, Osborne, J.P., et al. 2005, GCN Circ. 

3191 

Zeh, A., Klose, S., & Kann, D.A. 2006, ApJ, 637, 889 
Zhang, B., & Meszaros, P. 2001, ApJ, 552, L35 
Zhang, B., & Meszaros, P. 2002, ApJ, 571, 876 
Zhang, B., Fan, Y.Z., Dyks, J., et al. 2006, ApJ, 642 ,354 



